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Phase coexistence of gradient Gibbs states 


Abstract. We consider the (scalar) gradient fields p = (pb )—with b denoting the nearest- 
neighbor edges in Z 2 —that are distributed according to the Gibbs measure proportional 
to e~ / 3 H ^u(dp). Here H = V(pb) is the Hamiltonian, V is a symmetric potential, 
P > 0 is the inverse temperature, and v is the Lebesgue measure on the linear space defined 
by imposing the loop condition pb 1 + Pb 2 — + t?f> 4 for eac h plaquette (b 1, &2,63,64) 

in Z 2 . For convex V , Funaki and Spohn have shown that ergodic infinite-volume Gibbs 
measures are characterized by their tilt. We describe a mechanism by which the gradient 
Gibbs measures with non-convex V undergo a structural, order-disorder phase transition at 
some intermediate value of inverse temperature p. At the transition point, there are at least 
two distinct gradient measures with zero tilt, i.e., Epb = 0. 


1. Introduction 

1.1. Gradient fields 

One of the mathematical challenges encountered in the study of systems exhibit¬ 
ing phase coexistence is an efficient description of microscopic phase boundaries. 
Here various levels of detail are in general possible: The finest level is typically as¬ 
sociated with a statistical-mechanical model (e.g., a lattice gas) in which both the 
interface and the surrounding phases are represented microscopically; at the coars¬ 
est level the interface is viewed as a macroscopic (geometrical) surface between 
two structureless bulk phases. An intermediate approach is based on effective 
(and, often, solid-on-solid) models, in which the interface is still microscopic— 
represented by a stochastic field—while the structural details of the bulk phases 
are neglected. 

A simple example of such an effective model is a gradient field. To define this 
system, we consider a finite subset A of the d-dimensional hypercubic lattice 1r 
and, at each site of A and its external boundary <9/1, we consider the real-valued 
variable (f> x representing the height of the interface at x. The Hamiltonian is then 
given by 
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where the sum is over unordered nearest-neighbor pairs (. x , y). A standard example 
is the quadratic potential V(rf) = - nr/ 2 with k > 0; in general V is assumed to 
be a smooth, even function with a sufficient (say, quadratic) growth at infinity. The 
Gibbs measure takes the usual form 


P A (d<j)) = Z~ X Q- pHA ^d(j), (1.2) 

where d <fi is the |/1|-dimensional Lebesgue measure (the boundary values of <j> 
remain fixed and implicit in the notation), (3 > 0 is the inverse temperature and Z 
is a normalization constant. 

A natural question to ask is what are the possible limits of the Gibbs measures 
P/i(d</>) as A f Z d . Unfortunately, in dimensions d = 1,2, the fields (4> x ) x qa are 
very “rough” no matter how tempered the boundary conditions are assumed to be. 
As a consequence, the family of measures ( Pa) acj/' is not tight and no meaning¬ 
ful object is obtained by taking the limit A | Z d —i.e., the interface is delocalized. 
On the other hand, in dimensions d > 3 the fields are sufficiently smooth to permit 
a non-trivial thermodynamic limit—the interface is localized. These facts are es¬ 
tablished by combinations of Brascamp-Lieb inequality techniques and/or random 
walk representation (see, e.g., [16]) which, unfortunately, apply only for convex 
potentials with uniformly positive curvature. Thus, somewhat surprisingly, even 
for V(ij) = r / 4 the problem of localization in high-dimension is still open [24, 
Open Problem 1], 

As it turns out, the thermodynamic limit of the measures Pa is significantly 
less singular once we restrict attention to the gradient variables rj = (pb). These 
are defined by rjb = 4> y — <j> x where b is the nearest-neighbor edge (x, y) oriented 
in one of the positive lattice directions. Indeed, the 77 -marginal of P A (d<j>) always 
has at least one (weak) limit “point” as A —> Z . The limit measures satisfy a 
natural DLR condition and are therefore called gradient Gibbs measures. (Precise 
definitions will be stated below or can be found in [16,23].) One non-standard 
aspect of the gradient variables is that they have to obey a host of constraints. 
Namely, 

VG + r]b 2 = Vb 3 + tlbi ( 1 - 3 ) 

holds for each lattice plaquette (■ b \, 62 , 63 , 64 ), where the edges bj are listed coun¬ 
terclockwise and are assumed to be positively oriented. These constraints will be 
implemented at the level of a priori measure, see Sect. 2. 

It would be natural to expect that the character (and number) of gradient Gibbs 
measures depends sensitively on the potential V. However, this is not the case for 
the class of uniformly strictly-convex potentials (i.e., the V’s such that V"(p) > 
c_ > 0 for all 77 ). Indeed, Funaki and Spohn [17] showed that, in these cases, the 
translation-invariant, ergodic, gradient Gibbs measures are completely character¬ 
ized by the tilt of the underlying interface. Here the tilt is a vector u £ M. d such 
that 

Er]b = u ■ b (1.4) 

for every edge b —which we regard as a vector in R d . Furthermore, the correspon¬ 
dence is one-to-one, i.e., for each tilt there exists precisely one gradient Gibbs 
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measure with this tilt. Alternative proofs permitting extensions to discrete gradient 
models have appeared in Sheffield’s thesis [23]. 

It is natural to expect that a serious violation of the strict-convexity assump¬ 
tion on V may invalidate the above results. Actually, an example of a gradient 
model with multiple gradient Gibbs states of the same tilt has recently been pre¬ 
sented [23]; unfortunately, the example is not of the type considered above be¬ 
cause of the lack of translation invariance and its reliance on the discreteness of 
the fields. The goal of this paper is to point out a general mechanism by which 
the model ( 1 . 1 ) with a sufficiently non-convex potential V fails the conclusions of 
Funaki-Spohn’s theorems. 

1.2. Potentials of interest 

The mechanism driving our example will be the occurrence of a structural sur¬ 
face phase transition. To motivate the forthcoming considerations, let us recall 
that phase transitions typically arise via one of two mechanisms: either due to the 
breakdown of an internal symmetry, or via an abrupt turnover between energeti¬ 
cally and entropically favored states. The standard examples of systems with these 
kinds of phase transitions are the Ising model and the 7 -state Potts model with a 
sufficiently large q , respectively. In the former, at sufficiently low temperatures, 
there is a spontaneous breaking of the symmetry between the plus and minus spin 
states; in the latter, there is a first-order transition at intermediate temperatures be¬ 
tween q ordered, low-temperature states and a disordered, high-temperature state. 

Our goal is to come up with a potential V that would mimic one of the above 
situations. In the present context an analogue of the Ising model appears to be a 
double-well potential of the form, e.g., 

V(v) = ~? 7 *) 2 . ( 1 . 5 ) 

Unfortunately, due to the underlying plaquette constraints (1.3), the symmetry be¬ 
tween the wells cannot be completely broken and, even at the level of ground 
states, the system appears to be disordered. On 1? this can be demonstrated ex¬ 
plicitly by making a link to the ice model , which is a special case of the six vertex 
model [1]. A similar equivalence has been used [2] to study a roughening transition 
in an SOS interface. 

To see how the equivalence works exactly, note that the ground states of the 
system (1.5) are such that all 77 ’s equal ± 77 *. Let us associate a unit flow with each 
dual bond whose sign is determined by the value of i]b for its direct counterpart b. 
The plaquette constraint (1.3) then translates into a no-source-no-sink condition 
for this flow. If we mark the flow by arrows, the dual bonds at each plaquette are 
constrained to one of six zero-flux arrangements of the six vertex model; cf Fig. 1 
and its caption. The weights of all zero-flux arrangements are equal; we thus have 
the special case corresponding to the ice model. The ice model can be “exactly 
solved” [1]: The ground states have a non-vanishing residual entropy [22] and are 
disordered with infinite correlation length [1, Sect. 8 .10.III]. However, it is not 
clear how much of this picture survives to positive temperatures. 
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Fig. 1. The six plaquette configurations of minimal energy for the potential (1.5) on Z 2 and 
their equivalent ice model configurations at the corresponding vertex on the dual grid. The 
sign marks represent the signs of r/b along the side of the plaquette ( bi , b 2, &3,64)— with 
horizontal bonds 61, 63 oriented to the right and vertical bonds 62, &4 oriented upwards. The 
unit flow represented by the arrows runs upwards (downwards) through horizontal bonds 
with positive (negative) sign, and to the left (right) through vertical bonds with positive 
(negative) sign. The loop condition (1.3) makes the flow conserved (i.e., no sources or 
sinks). 

The previous discussion shows that it will be probably quite hard to realize 
a symmetry-breaking transition in the context of the gradient model (1.1). It is 
the order-disorder mechanism for phase transitions that seems considerably more 
promising. There are two canonical examples of interest: a potential with two cen¬ 
tered wells and a triple-well potential ; see Fig. 2. Both of these lead to a gradient 
model which features a phase transition, at some intermediate temperature, from 
states with the 77 ’s lying (mostly) within the thinner well to states whose 77 ’s fluc¬ 
tuate on the scale of the thicker well(s). 

Our techniques apply equally to these—as well as other similar—cases pro¬ 
vided the widths of the wells are sufficiently distinct. Notwithstanding, the analysis 
becomes significantly cleaner if we abandon temperature as our principal parame¬ 
ter (e.g., we set f3 = 1) and consider potentials V that are simply defined by 

e -T(0 = p Q — Korj 2 /2 + (1 _ p ) e -/W/2_ (1.6) 

Here n 0 and k d are positive numbers and p is a parameter taking values in [0,1], 
For appropriate values of the constants, V defined this way will have a graph as in 
Fig. 2(a). To get the graph in part (b), we would need to consider V’s of the form 

e ~ v (y) = ps~ K °y 2 / 2 + ^ ~ p e ~ K v(y-v*) 2 / 2 _|_ ^ ~ p e -K D (v+y*) 2 / 2 ^ (1.7) 
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Fig. 2. Two canonical examples of potentials that will lead to a structural surface phase 
transition. The picture labeled (a) is obtained by superimposing—in the sense of (1.6)— 
two symmetric wells of (significantly) different widths. Part (b) of the figure represents the 
triple-well potential as defined in (1.7). For the application of our technique of proof, it only 
matters that the widths of the wells are sufficiently different. 


where ±/y* are the (approximate) locations of the off-center wells. 

The idea underlying the expressions (1.6) and (1.7) is similar to that of the 
Fortuin-Kasteleyn representation of the Potts model [12]. In the context of conti¬ 
nuous-spin models similar to ours, such representation has fruitfully been used by 
Zahradnfk [25]. Focusing on (1.6), we can interpret the terms on the right-hand side 
of (1.6) as two distinct states of each bond. (We will soon exploit this interpretation 
in detail.) The indexing of the coupling constants suggests the names: “O” for 
ordered and “D” for disordered. 

It is clear that the extreme values of p (near zero or near one) will be domi¬ 
nated by one type of bonds; what we intend to show is that, for n 0 and k d suf¬ 
ficiently distinct from each other, the transition between the “ordered" and “dis¬ 
ordered" phases is (strongly) first order. Similar conclusions and proofs—albeit 
more complicated—apply also to the potential (1.7). However, for clarity of expo¬ 
sition, we will focus on the potential (1.6) for the rest of the paper (see, however. 
Sect. 2.5). In addition, we will also restrict ourselves to two dimensions, even 
though the majority of our results are valid for all d> 2. 


2. Main results 

2.1. Gradient Gibbs measures 


We commence with a precise definition of our model. Most of the work in this 
paper will be confined to the lattice torus of L x L sites in Z 2 , so we will start 
with this particular geometry. Choosing the natural positive direction for each lat¬ 
tice axis, let B/, denote the corresponding set of positively oriented edges in T /,. 
Given a configuration {<j) x )xeT L , we introduce the gradient field rj = V</> by as¬ 
signing the variable = <j> y — <f> x to each b = ( x , y) g 1^. The product Lebesgue 
measure Ylxjto induces a (er-finite) measure on the space R Bi via 





( 2 . 1 ) 
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where S denotes the Dirac point-mass at zero. 

We interpret the measure vl as an a priori measure on gradient configura¬ 
tions 77 € R Bi . Since the 77’s arise as the gradients of the ( p’s it is easy to check 
that vl is entirely supported on the linear subspace Xl C R Bl of configurations 
determined by the condition that the sum of signed 77’s—with a positive or neg¬ 
ative sign depending on whether the edge is traversed in the positive or negative 
direction, respectively—vanishes around each closed circuit on 1 ’ L . (Note that, in 
addition to ( 1 . 3 ), the condition includes also loops that wrap around the torus.) We 
will refer to such configurations as curl-free. 

Next we will define gradient Gibbs measures on T/. For later convenience 
we will proceed in some more generality than presently needed: Let {Vb)b^m L be 
a collection of measurable functions 14 : R —> [ 0 , 00 ) and consider the partition 
function 


Z L,(V b ) 



exp 


{- Vb (Vb)}vL{dr)). 
bea L 


( 2 . 2 ) 


Clearly, .Zl.( t4) > 0 and, under the condition that 77 1—> e Vh( - r D is integrable 
with respect to the Lebesgue measure on R, also Z L iy b \ < 00 . We may then 
define PL,(v b ) to be the probability measure on RI Bl I given by 


p l,(\ 4 )( d? ?) 


J L,(v b ) 


ex p{- v b(Vb)} t> L {drf). 


b 


( 2 . 3 ) 


This is the gradient Gibbs measure on Tl corresponding to the potentials ( 14 ). In 
the situations when 14 = V for all b —which is the principal case of interest in this 
paper—we will denote the corresponding gradient Gibbs measure on Tl by Pl,v- 
It is not surprising that PL,(v b ) obeys appropriate DLR equations with respect 
to all connected A C Tl containing no topologically non-trivial circuit. Explic¬ 
itly, if 77,4c in A c is a curl-free boundary condition, then the conditional law of 77 4 
given 77,4c is 


PL,(v b ){driA\VA*) 


Za{jia<) 


ex p{- ^2 H(%)} VAidriAlVA*)- 


beA 


( 2 . 4 ) 


Here PL,(v b ){dVA\VA°) is the conditional probability with respect to the (tail) er- 
algebra . 4 /i generated by the fields on A c , Za{t]a c ) is the partition function in A, 
and t'A (d?7yi | tja c ) is the a priori measure induced by vl on t/a given the boundary 
condition 77 a c ■ 

As usual, this property remains valid even in the thermodynamic limit. We 
thus say that a measure on p is an infinite-volume gradient Gibbs measure if it 
satisfies the DLR equations with respect to the specification ( 2 . 4 ) in any finite 
set 4 C Z 2 . (As is easy to check—e.g., by reinterpreting the 77’s back in terms 
of the </>’s—77/1 (d?7/t 1 77,4c) is independent of the values of ?;> outside any circuit 
winding around A, and so it is immaterial that it originated from a measure on 
torus.) 
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An important aspect of our derivations will be the fact that our potential V 
takes the specific form ( 1 . 6 ), which can be concisely written as 


e -v(v) 


g(d K)e~? Kri \ 


(2.5) 


where g is the probability measure g = p6 Ko + (1 — p)S KD . It follows that the 
Gibbs measure Pl.v can be regarded as the projection of the extended gradient 
Gibbs measure. 


Qt(dt?)dft) = ex p{-^ KbT lb} VL(dv)eL(dK), (2.6) 


be M l 


to the cr-algebra generated by the ry’s. Here g^ is the product of measures g, one 
for each bond in Ml. As is easy to check, conditioning on (pb, Kb)beA c yields 
the corresponding extension QA{dpAdKA\r]A c ) of the finite-volume specification 
(2.4)—the result is independent of the k’s outside A because, once p A c is fixed, 
these have no effect on the configurations in A. 

The main point of introducing the extended measure is that, if conditioned on 
the k’s, the variables pb are distributed as gradients of a Gaussian field—albeit with 
a non-translation invariant covariance matrix. As we will see, the phase transition 
proved in this paper is manifested by a jump-discontinuity in the density of bonds 
with Kb = Ko which at the level of 77 -marginal results in a jump in the characteristic 
scale of the fluctuations. 


Remark 2.1. Notably, the extended measure Ql plays the same role for Pl.v as 
the so called Edwards-Sokal coupling measure [10] does for the Potts model. Sim¬ 
ilarly as for the Edwards-Sokal measures [3,18], there is a one-to-one correspon¬ 
dence between the infinite-volume measures on p's and the corresponding infinite- 
volume extended gradient Gibbs measures on (p, k)’s. Explicitly, if // is an infinite- 
volume gradient Gibbs measure for potential V, then fi, defined by (extending the 
consistent family of measures of the form) 

£((%, Kb)b£A & Ax B) = [ g A (dK) eJa a T\ e-l'W+n^ A > (2 .7) 

J » V b£A > 

is a Gibbs measure with respect to the extended specifications Qa{’\va c )- For the 
situations with only a few distinct values of Kb, it may be of independent interest 
to study the properties of the rc-marginal of the extended measure, e.g., using the 
techniques of percolation theory. However, apart from some remarks in Sect. 2.3, 
we will not pursue these matters in the present paper. 


2.2. Phase coexistence of gradient measures 

Now we are ready to state our main results. Throughout we will consider the po¬ 
tentials V of the form (1.6) with k 0 ^ Kd- As a moment’s thought reveals, the 
model is invariant under the transformation 
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for any fixed 9 0. In particular, without loss of generality, one could assume 

from the beginning that KoKd = 1 and regard k °/ Kd as the sole parameter of the 
model. However, we prefer to treat the two terms in (1.6) on an equal footing, and 
so we will keep the coupling strengths independent. 

Given a shift-ergodic gradient Gibbs measure, recall that its tilt is the vector u 
such that (1.4) holds for each bond. The principal result of the present paper is the 
following theorem: 

Theorem 2.2. For each e > 0 there exists a constant c = c(e) >0 and, if 


K 0 > ck d , (2.9) 

a number p t € (0,1) such that, for interaction V with p = p t , there are two 
distinct, infinite-volume, shift-ergodic gradient Gibbs measures /i or d and pais of 
zero tilt for which 

fiord(|??fc| > -2=) < e+ VA > 0, (2.10) 

and 

Mdis(W < -y==) - e + ciA 1/4 , VA > 0. (2.11) 

Here Ci is a constant of order unity. 


Remark 2.3. An inspection of the proof actually reveals that the above bounds are 
valid for any e satisfying e > C 2 ( Kd / Ko ) 1//8 , where Ci is a constant of order unity. 


As already alluded to, this result is a consequence of the fact that the density 
of ordered bonds, i.e., those with Kb = Ko, undergoes a jump at p = p t . On the 
torus, we can make the following asymptotic statements: 

Theorem 2.4. Let R°f d denote the fraction of ordered bonds on TV, i.e., 

^ rd =MlEW«o}- (2-12) 

1 L 1 be M l 

For each e > 0 there exists c = c(e) > 0 such that the following holds: Under the 
condition (2.9), and for p t as in Theorem 2.2, 

lim Q L (R°f d <e) = l, p < p t (2.13) 

L— >oo 

and 

lim Q L (R°l A > 1 — e) = 1, p > p t . (2.14) 

L —►oo 

The present setting actually permits us to determine the value of p t via a du¬ 
ality argument. This is the only result in this paper which is intrinsically two- 
dimensional (and intrinsically tied to the form (1.6) of V ). All other conclusions 
can be extended to d > 2 and to more general potentials. 


Theorem 2.5. Let d = 2. If k °/ Kd 1, then p t is given by 


Pt 

1 - Pt 



(2.15) 


Theorem 2.4 is proved in Sect. 4.2, Theorem 2.2 is proved in Sect. 4.3 and 
Theorem 2.5 is proved in Sect. 5.3. 
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2.3. Discussion 

The phase transition described in the above theorems can be interpreted in several 
ways. First, in terms of the extended gradient Gibbs measures on toms, it clearly 
corresponds to a transition between a state with nearly all bonds ordered (Kb = Ko) 
to a state with nearly all bonds disordered (Kb = k d ). Second, looking back at the 
inequalities ( 2 . 10 - 2 . 11 ), most of the if s will be of order at most in the 

ordered state while most of them will be of order at least l/ifitH in the disor¬ 
dered state. Hence, the corresponding (effective) interface is significantly rougher 
at p < p { than it is at p > p t (both phases are rough according to the standard 
definition of this term) and we may thus interpret the above as a kind of first-order 
roughening transition that the interface undergoes at p t . Finally, since the gradient 
fields in the two states fluctuate on different characteristic scales, the entropy (and 
hence the energy) associated with these states is different; we can thus view this 
as a standard energy-entropy transition. (By the energy we mean the expectation 
of V ( 77 b); notably, the expectation of KbP^ is the same in both measures; cf (4.35).) 
Energy-entropy transitions for spin models have been studied in [9,20,21] and, 
quite recently, in [ 11 ], 

Next let us turn our attention to the conclusions of Theorem 2.4. We actually 
believe that the dichotomy (2.13-2.14) applies (in the sense of almost-sure limit 
of R°l A as L —> oo) to all translation-invariant extended gradient Gibbs states with 
zero tilt. The reason is that, conditional on the k’s, the gradient fields are Gaus¬ 
sian with uniformly positive stiffness. We rest assured that the techniques of [17] 
and [23] can be used to prove that the gradient Gibbs measure with zero tilt is 
unique for almost every configuration of the k’s; so the only reason for multi¬ 
plicity of gradient Gibbs measures with zero tilt is a phase transition in the k- 
marginal. However, a detailed write-up of this argument would require developing 
the precise—and somewhat subtle—correspondence between the gradient Gibbs 
measures of a given tilt and the minimizers of the Gibbs variational principle 
(which we have, in full detail, only for convex periodic potentials [23]). Thus, 
to keep the paper at manageable length, we limit ourselves to a weaker result. 

The fact that the transition occurs at p t satisfying (2.15) is a consequence of a 
duality between the k- marginals at p and 1 — p. More generally, the duality links 
the marginal law of the configuration (Kb) with the law of (1 /Kb)\ see Theorem 5.3 
and Remark 5.4. [At the level of gradient fields, the duality provides only a vague 
link between the flow of the weighted gradients ( v / ~KbPb ) along a given curve and 
its flux through this curve. Unfortunately, this link does not seem to be particularly 
useful.] The point p - p t is self-dual which makes it the most natural candidate 
for a transition point. It is interesting to ponder about what happens when k °/ Kd 
decreases to one. Presumably, the first-order transition (for states at zero tilt) dis¬ 
appears before k °/ Kd reaches one and is replaced by some sort of critical behavior. 
Here the first problem to tackle is to establish the absence of first-order phase tran¬ 
sition for small k °/ Kd — 1. Via a standard duality argument (see [8]) this would 
yield a power-law lower bound for bond connectivities at p t . 

Another interesting problem is to determine what happens with measures of 
non-zero tilt. We expect that, at least for moderate values of the tilt u , the first- 
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order transition persists but shifts to lower values of p. Thus, one could envision 
a whole phase diagram in the p-u plane. Unfortunately, we are unable to make 
any statements of this kind because the standard ways to induce a tilt on the torus 
(cf [17]) lead to measures that are not reflection positive. 


2.4. Outline of the proof 


We proceed by an outline of the principal steps of the proof to which the remainder 
of this paper is devoted. The arguments are close in spirit to those in [9,20,21]; 
the differences arise from the subtleties in the setup due to the gradient nature of 
the fields. 

The main line of reasoning is basically thermodynamical: Consider the k- 
marginal of the extended torus state Ql which we will regard as a measure on 
configurations of ordered and disordered bonds. Let \{p) denote (the L —» oo 
limit of) the expected fraction of ordered bonds in the torus state at parameter p. 
Clearly x(p) increases from zero to one as p sweeps through [0,1]. The princi¬ 
pal observation is that, under the assumption k °/ Kd 1, the quantity \(1 — %) is 
small, uniformly in p. Hence, p i—> x(p) must undergo a jump from values near 
zero to values near one at some p t £ (0,1). By usual weak-limiting arguments 
we construct two distinct gradient Gibbs measures at p t , one with high density of 
ordered bonds and the other with high density of disordered bonds. 

The crux of the matter is thus to justify the uniform smallness of \(1 —x). This 
will be a consequence of the fact that the simultaneous occurrence of ordered and 
disordered bonds at any two given locations is (uniformly) unlikely. For instance, 
let us estimate the probability that a particular plaquette has two ordered bonds 
emanating out of one corner and two disordered bonds emanating out of the other. 
Here the technique of chessboard estimates [13-15] allows us to disseminate this 
pattern all over the torus via successive reflections (cf Theorem 4.2 in Sect. 4.1). 
This bounds the quantity of interest by the 1 //7-power of the probability that 
every other horizontal (and vertical) line is entirely ordered and the remaining 
lines are disordered. The resulting “spin-wave calculation”—i.e., diagonalization 
of a period-2 covariance matrix in the Fourier basis and taking its determinant—is 
performed (for all needed patterns) in Sect. 3. 

Once the occurrence of a “bad pattern” is estimated by means of various spin- 
wave free energies, we need to prove that these “bad-pattern” spin-wave free ener¬ 
gies are always worse off than those of the homogeneous patterns (i.e., all ordered 
or all disordered)—this is the content of Theorem 3.3. Then we run a standard 
Peierls’ contour estimate whereby the smallness of x(l — x) follows. Extracting 
two distinct, infinite-volume, ergodic gradient Gibbs states p , on j and pa s at p = p t , 
it remains to show that these are both of zero tilt. Here we use the fact that, condi¬ 
tional on the ft’s, the torus measure is symmetric Gaussian with uniformly positive 
stiffness. Hence, we can use standard Gaussian inequalities to show exponential 
tightness of the tilt, uniformly in the ft’s; cf Lemma 4.8. Duality calculations (see 
Sect. 5) then yield p = p t . 
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2.5. Generalizations 

Our proof of phase coexistence applies to any potential of the form shown in 
Fig. 2—even if we return to parametrization by 3. The difference with respect 
to the present setup is that in the general case we would have to approximate the 
potentials by a quadratic well at each local minimum and, before performing the 
requisite Gaussian calculations, estimate the resulting errors. 

Here is a sketch of the main ideas: We fix a scale A and regard pb to be in 
a well if it is within A of the corresponding local minimum. Then the requisite 
quadratic approximation of /3-times energy is good up to errors of order 3A 3 . The 
rest of the potential “landscape” lies at energies of at least order A 2 and so it will 
be only “rarely visited” by the 77 ’s provided that 3A 2 1. On the other hand, the 
same condition ensures that the spin-wave integrals are essentially not influenced 
by the restriction that rjb be within A of the local minimum. Thus, to make all 
approximations work we need that 

3 A 3 < 1 < 3 A 2 (2.16) 

which is achieved for /3 1 by, e.g., A = 3~^ ■ This approach has recently been 

used to prove phase transitions in classical [4,5] as well as quantum [ 6 ] systems 
with highly degenerate ground states. We refer the reader to these references for 
further details. 

A somewhat more delicate issue is the proof that both coexisting states are of 
zero tilt. Here the existing techniques require that we have some sort of uniform 
convexity. This more or less forces us to use the V’s of the form 

V(rj) = - logf^e-^), (2.17) 

'3 ' 

where the Vfs are uniformly convex functions. Clearly, our choice (1.6) is the 
simplest potential of this type; the question is how general the potentials obtained 
this way can be. We hope to return to this question in a future publication. 

3. Spin-wave calculations 

As was just mentioned, the core of our proofs are estimates of the spin-wave free 
energy for various regular patterns of ordered and disordered bonds on the torus. 
These estimates are rather technical and so we prefer to clear them out of the way 
before we get to the main line of the proof. The readers wishing to follow the 
proof in linear order may consider skipping this section and returning to it only 
while reading the arguments in Sect. 4.2. Throughout this and the forthcoming 
sections we assume that L is an even integer. 

3.1. Constrained partition functions 

We will consider six partition functions Z 10 , Zl, iren Zl,ud, Zl.mp and 

Zl.ma °n T l that correspond to six regular configurations each of which is ob¬ 
tained by reflecting one of six possible arrangements of “ordered” and “disor¬ 
dered” bonds around a lattice plaquette to the entire torus. These quantities will be 
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Fig. 3. Six possible arrangements of “ordered” and “disordered” bonds around a lattice 
plaquette. Here the “ordered” bonds are represented by solid lines and the “disordered” 
bonds by wavy lines. Each inhomogeneous pattern admits other rotations which are not 
depicted. The acronyms stand for: (top row) Ordered. Disordered and U-shape Ordered and 
(bottom row) U-shape Disordered. Mixed Periodic and Mixed Aperiodic, respectively. 


the “building blocks” of our analysis in Sect. 4. The six plaquette configurations 
are depicted in Fig. 3. 

We begin by considering the homogeneous configurations. Here Zl,o is the 
partition function Z L i Vb , for all edges of the “ordered” type: 

V b {r]) = -logp+ be M l . (3.1) 

Similarly, Zi n is the quantity Zj Jt (v b ) for 

Vb(rf) = - log(l ~P) + b e M l , (3.2) 

i.e., with all edges “disordered.” 

Next we will define the partition functions Z/uo and Zl,vd which are ob¬ 
tained by reflecting a plaquette with three bonds of one type and the remaining 
bond of the other type. Let us split B l into the even B^ en and odd B^ dd horizontal 
and vertical edges—with the even edges on the lines of sites in the x direction 
with even y coordinates and lines of sites in y direction with even x coordinates. 
Similarly, we will also consider the decomposition of B^ into the set of horizontal 
edges B^ or and vertical edges B^ rt . Letting 


Vb{v) 


log p+ \k 0 r] 2 , if be B dor U B e L ven 

log(l — p) + jKd?) 2 , otherwise, 


the partition function Zl, uo then corresponds to the quantity Z^ ( Vhr The partition 
function Zl,\j d is obtained similarly; with the roles of “ordered” and “disordered” 
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interchanged. Note that, since we are working on a square toms, the orientation of 
the pattern we choose does not matter. 

It remains to define the partition functions Z/yp and Z l,ma corresponding to 
the patterns with two “ordered” and two “disordered” bonds. For the former, we 
simply take Zj J ^y b ' ) with the potential 


Mv) 


log p+±K 0 r] 2 , ifbem h ° r , 

log(l -p) + ±kW, if 6 G B™ rt . 


Note that the two types of bonds are arranged in a “mixed periodic” pattern; hence 
the index MP. As to the quantity Z^ ma. here we will consider a “mixed aperiodic” 
pattern. Explicitly, we define 


Vb(v) 


-logp + \norj 2 , 
log(l p) + \k-dV 2 i 


if be B™, 

if b G l^ dd . 


(3.5) 


The “mixed aperiodic” partition function Zl ma is the quantity Z L j Vh) for this 
choice of (14). Again, on a square torus it is immaterial for the values of Zl,mp 
and which orientation of the initial plaquette we start with. 

As usual, associated with these partition functions are the corresponding free 
energies. In finite volume, these quantities can be defined in all cases by the for¬ 
mula 

FlAp ) = -^2 log ( 2 ^%°-t) » a = 0 ,D,UO,UD,MP,MA, (3.6) 

where the factor (27r)tK L ^ has been added for later convenience and where 
the p-dcpcndence arises via the corresponding formulas for 14 in each particular 
case. 


3.2. Limiting free energies 

The goal of this section is to compute the thermodynamic limit of the Fl ]Q ’s. For 
homogeneous and isotropic configurations, an important role will be played by the 
momentum representation of the lattice Laplacian D(k) = |1 — e lfel | 2 +11 — e lk2 1 2 
defined for all k = (k \, Ap) in the corresponding Brillouin zone k £ [— n, tt\ x 
[—7r, 7r]. Using this quantity, the “ordered" free energy will be simply 

II' dk 

Fo(p) = 2 log p — I —2 log {n 0 D{k)}, (3.7) 

1 J[-n,7T] 2 (Z7T) 

while the disordered free energy boils down to 

1 r fife ^ 

F d (p) = -21og(l -p) + - / —— ^logjACo D{k)}. (3.8) 

2 J [—7T,7r] 2 (27 t)“ 
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It is easy to check that, despite the logarithmic singularity at k = 0, both integrals 
converge. The bond pattern underlying the quantity Z/ M p lacks rotation invari¬ 
ance and so a different propagator appears inside the momentum integral: 


Fmp(p) = -log[p(l-p)] 

r\U 

p log{K 0 |l - e lfcl | 2 + k d |1 - e lfc2 | 2 }. (3.9) 

Again, the integral converges as long as (at least) one of k 0 and re D is strictly 
positive. 

The remaining partition functions come from configurations that lack transla¬ 
tion invariance and are “only” periodic with period two. Consequently, the Fourier 
transform of the corresponding propagator is only block diagonal, with two or four 
different fc’s “mixed” inside each block. In the UO cases we will get the function 

1 i r dfc 

^uo(p) = --log[p 3 (l -p)] + - / —^ log{det7T uo (fc)}, (3.10) 

Z 4 J[- 7T,7r] 2 (A 7r ) 

where II\jo(k) is the 2 x 2-matrix 



^uo(fe) = 


«o|a -| 2 + b(no + kd)|&-| 


|2 _L i/V- -L I 2 i(«; 0 — k d )|6_| 2 


- k d )|6_| 


^o | a +| 2 + Kv)\b- 


with a± and b± defined by 

a± = l±e ifcl and &± = l±e ifc2 . 


(3.11) 

(3.12) 


The extra factor V 2 —on top of the usual V 2 —in front of the integral arises because 
det TZuo (fe) combines the contributions of two Fourier models; namely k and k + 
7rei. A calculation shows 

detTTuo(fc) > K 0 2 |o_| 2 |a+| 2 + k 0 «T>IM 4 > (3.13) 

implying that the integral in (3.10) converges. The free energy Fud is obtained by 
interchanging the roles of no and kd and of p and (1 — p). 

In the MA-cases we will assume that kq ^ kd— otherwise there is no dis¬ 
tinction between any of the six patterns. The corresponding free energy is then 
given by 


Fma(p) 


log[p(l ~p)\ 




d k 

w? 


log| ( h ° 2 K ° ) det77 M A(fc)|. 


Here IlMA{k) is the 4 x 4-matrix 


(3.14) 


n ma(^) 


fr(\a_\ 2 + \b.\ 2 ) \b-\ 2 |«-| 2 0 \ 

|6_| 2 r(|a + | 2 + |6_| 2 ) 0 |a + | 2 

-I 2 0 r(|a_| 2 + |6 + | 2 ) \b + \ 2 

0 |a + | 2 \b+\ 2 r(|a + | 2 + |6+| 2 )/ 

(3.15) 


V 
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with the abbreviation 


Ko + k d 

r = -. 

«o - «D 


(3.16) 


Note that r > 1 in the cases of our interest. 

Observe that detflMA(fc) is a quadratic polynomial in r 2 , i.e., detTTMA = 
Ar 4 + Br 2 + C. Moreover, IIuAik) annihilates (1, —1, —1,1) when r = 1, and 
so r 2 = 1 is a root of Ar 4 + Br 2 + C. Hence det IlMA(k) = (r 2 — 1 ){Ar 2 — C), 
i.e.. 


det77 M A(fc) = (r 2 - 1)j-(|a+| 2 |a_| 2 - |6 + | 2 |fe_| 2 ) 2 
+ (|a+| 2 + |6 + | 2 )(|a_| 2 + |fc+| 2 )(|a+| 2 + |6_| 2 )(|«_| 2 + |6_| 2 )r 2 }. 

(3.17) 

Setting r = 1 inside the large braces yields 

detf7 M A(fc) > 4(r 2 - l)|a_| 2 |a+| 2 |b_| 2 |6+| 2 , (3.18) 

implying that the integral in (3.14) is well defined and finite. 

Remark 3.1. The fact that 11 w\(k) has zero eigenvalue at r = 1 is not surprising. 
Indeed, r = 1 corresponds to « D = 0 in which case a quarter of all sites in the MA- 
pattern get decoupled from the rest. This indicates that the partition function blows 
up (at least) as (r — l)~l Ti l/ 4 as r | 1 implying that there should be a zero 
eigenvalue at r = 1 per each 4 x 4-block TImaC^)- 

A formal connection between the quantities in (3.6) and those in (3.7-3.14) is 
guaranteed by the following result: 

Theorem 3.2. For all a = O, D, UO, UD, MP, MA and uniformly in p € (0,1), 

lim F Lta (p) = F a (p). (3.19) 

L —»oo 

Proof. This is a result of standard calculations of Gaussian integrals in momentum 
representation. We begin by noting that the Lebesgue measure d@ x can be 
regarded as the product of v r , acting only on the gradients of ()>, and d f z for some 
fixed z &Tl. Neglecting temporarily the a priori bond weights p and (1 — p), the 
partition function Zf. rt , a = O, D, UO, UD, MP, MA, is thus the integral of the 
Gaussian weight (27r|Ti|) _1//2 e~5( < ^> c 'Q F) against the measure JX r df x , where 
the covariance matrix C a is defined by the quadratic form 

(0,67-V) = ]T K ( b a \\/ b f) 2 + (3-20) 

6£B l I L ' x€Tl 

Here (k^) are the bond weights of pattern a. Indeed, the integral over d <j> z with 
the gradient variables fixed yields l) 1 ^ 2 which cancels the term in front of 

the Gaussian weight. The purpose of the above rewrite was to reinsert the “zero 
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mode” <po = |TV| -1 / 2 <p x into the partition function; <f >o was not subject to 
integration due to the restriction to gradient variables. 

To compute the Gaussian integral, we need to diagonalize C a . For that we will 
pass to the Fourier components <pk = |TV \~ 1 / 2 <$>xt lk x with the result 

(0,C-V)= E Mi’(s k ,o6v, 0 + E 

k,k'e T l <7=1,2 

(3.21) 

where TV = {x( n i> n2 ) : 0 < ni,ri 2 < L} is the reciprocal torus, 5 p q is the 
Kronecker delta and 


A 


(<V 

k,k' 


1 

irii 


E 


(a) 

K (x,x+e„) 


gi (k'-k)-x 


x£T L 


(3.22) 


Now if the horizontal part of (Kj '/’) is translation invariant in the 7 -th direction, 
then A^ k , = 0 whenever /c 7 V while if it is “only” 2 -periodic, then A^ k , = 0 
unless k 7 = kj or fc 7 = kj + 7rmod27r. Similar statements apply to the ver¬ 
tical part of («E) and A^ L. Since all of our partition functions come from 2- 
periodic configurations, the covariance matrix can be cast into a block-diagonal 
form, with 4x4 blocks O a (k) collecting all matrix elements that involve the 
momenta (fe, k + 7 rei, k + 7 re 2 , k + 7 rei + 7 re 2 ). Due to the reinsertion of the 
“zero mode”—cf (3.20)—all of these blocks are non-singular (see also the explicit 
calculations below). 

lienee we get that, for all a = O, D, UO, UD, MP, MA, 




- 


( 27r )i(L 2 -i) L p 


0--p) 


n d 


n 

fee Ti 


det 0 a {k) \ 


(3.23) 


where No and No denote the numbers of ordered and disordered bonds in the un¬ 
derlying bond configuration and where the exponent Vs takes care of the fact that 
in the product, each k gets involved in four distinct terms. Taking logarithms and 
dividing by |TV,|, the sum over the reciprocal torus converges to a Riemann inte¬ 
gral over the Brillouin zone [— it , 7t] x [— 7 r, 7r] (the integrand has only logarithmic 
singularities in all cases, which are harmless for this limit). 

It remains to justify the explicit form of the free energies in all cases under con¬ 
siderations. Here the situations a = O, D, MP are fairly standard, so we will focus 
on a = UO and a = MA for which some non-trivial calculations are needed. In 
the former case we get that 


A 


(i) 

k,k' 


Ko<Vc,fe'i 


A (2) _ tto + Kd 
A k,k ~ 0 


A (2) 

^ fc,fc+7rei 



(3.24) 


with A ( k a L = 0 for all values that are not of this type. Plugging into (3.21) we find 
that the (k, k + Tre^-subblock of 0\jo(k) reduces essentially to the 2 x 2 -matrix 
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in (3.11). Explicitly, 


0uo(fc) diag((ifc,o, ^ > <5fc,7re 2 > <5fc,-7rei+7re 2 ) d - 

Since whenever ^ 0, the block matrix Quo (k) will only be a 

function of moduli-squared of a± and 6_. Using (3.25) in (3.23) we get (3.10). 

As to the MA-case the only non-zero elements of are 


n\jo(k) 0 

0 n\jo(k 4 


7re 2 ) 
(3.25) 


A 1 ) _ tw __ 

ill, L - iiL L - 


( 2 ) 


tto + k D 


fc.fc 


fe.fc 


and 


A (1) = A (2) 

' ri fc,fe+7re2 ' ri fc,fe+7rei 



(3.26) 


So, again, k' a = k„ whenever A^ k , ^ 0 and so Qma(&) depends only on |a±| 2 
and \b±\ 2 . An explicit calculation shows that 


Qma(A;) = diag(<5 fc ,o,4 ,7rei ) Sk ,ne2 j Sk ,7rei+7re2 ) + ( K ° KD ) flMA(fc), (3.27) 
where TImaI^) is as in (3.15). Plugging into (3.23), we get (3.14). □ 


3.3. Optimal patterns 


Next we establish the crucial fact that the spin-wave free energies corresponding 
to inhomogeneous patterns UO, UD, MP, MA exceed the smaller of Fq and f j) by 
a quantity that is large, independent of p, once k 0 > Kd- 


Theorem 3.3. There exists q £ R such that if kd < £ no with £ G (0,1), then for 
all p G (0,1), 

min F a (p) - min F & (p) > - log — + - log(l - £) + ci. (3.28) 

a=UO,UD,MP,MA a=0,D 8 K D 4 


Proof. Let us use / and J to denote the integrals 

1 r dk 

I= n j log{D(fc)}, and J = 

2 J[-ir,nP l 2 ^) 



(3.29) 


We will prove (3.28) with c\ = J — I. 

First, we have 

Foip) = -2log p + - log k 0 + I 

and 

Fd(p) = -2 log(l -p) + - log k d + I, 
while an inspection of (3.14) yields 


(3.30) 

(3.31) 


F M a(p) = — log[p(l — p)\ 

If d k 

8 J[- W ,7r] 2 (2 tt)^ 


logtcottD(o - K D ) 2 |a+a_6 + 6_|“| 


> -logjp(l-p)] +^logK 0 + ilogK D + ilog(l-^) +J. (3.32) 
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Using that 

min{F 0 ,F D } <-(F 0 + F D ), (3.33) 

we thus get 

Fma(p) - min{F 0 (p), F D (p)} > i log — + i log(l - £) + J - I, (3.34) 

O hi D 4 

which agrees with (3.28) for our choice of ci. 

Coming to the free energy Fyo, using (3.13) we evaluate 

det7T uo (fc) > K 0 2 |a_| 2 |a+| 2 = (—) K 0 3/2 K V 1/2 \a- | 2 |o + | 2 (3.35) 

yielding 

Fuo(p) > ^ log[p 3 (1 ~p)\ + log — + ~log kq 4" ^logK D + J- (3.36) 

Z o Kd o o 

Bounding 

min{F 0 ,F D } < ^F 0 + ^F D (3.37) 

we thus get 

Fvoip) - min{F q {p ), F n (p)\ > ^ log — + J - I, (3.38) 

o Kq 

in agreement with (3.28). The computation for Fjjd is completely analogous, in¬ 
terchanging only the roles of k 0 and as well as p and (1 — p). From the lower 
bound 

det77 u0 (fc) > K 0 K D \b-\ 4 = (—') ' «o 1/2 Kd 3/2 |M 4 (3.39) 

and the inequality 

min{F 0 ,F D } < ^F 0 + jF d , (3.40) 

we get again 

Fvd(p) - min{F 0 (p), F D (p)} > ^ log — + J-I, (3.41) 

o Kq 

which is identical to (3.38). 

Finally, for the free energy Fmp, we first note that 

^MP(P) > - l0g[p(l -p)\+~ log Kq + J, (3.42) 

which yields 

F M p(p) - F D (p) > log -—- + i log — + J - I. (3.43) 

P 2 Kq 

Under the condition that log —^ > — | log we again get (3.28). For the com¬ 
plementary values of p, we will compare Fmp with Fo: 

Fvtp(p) - F 0 (p) > log —^-F J - I, (3.44) 

1 ~P 

Since we now have log —- < — | log this yields (3.28) with the above choice 
of Cl. 
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4. Proof of phase coexistence 

In this section we will apply the calculations from the previous section to the proof 
of Theorems 2.2 and 2.4. Throughout this section we assume that k 0 > «d anc * 
that L is even. We begin with a review of the technique of chessboard estimates 
which, for later convenience, we formulate directly in terms of extended configu¬ 
rations (Kb,T]b)- 


4.1. Review ofRP/CE technology 


Our principal tool will be chessboard estimates, based on reflection positivity. To 
define these concepts, let us consider the torus T p and let us split T p into two 
symmetric halves, Tp and T£, sharing a “plane of sites” on their boundary. We 
will refer to the set Tp fl Tp as plane of reflection and denote it by P. The half¬ 
tori Tp inherit the nearest-neighbor structure from Tl; we will use Bp to denote 
the corresponding sets of edges. On the extended configuration space, there is a 
canonical map Op: R® 1, x {kcKd}* 1, — »■ R® 1, x {koj^d}® 1, — induced by the 
reflection of Tp into Tp through P —which is defined as follows: If b,b' £ Mp 
are related via b' = Op(b), then we put 


(°pv) b = 


-Ob', 

Ob', 


if b J_ P, 

if Ml P, 


(4.1) 


and 

( OpK)b = Kb’■ (4.2) 

Here b L P denotes that b is orthogonal to p while b || P indicates that b is 
parallel to P. The minus sign in the case when b _L P is fairly natural if we recall 
that ?7b represents the difference of <f> x between the endpoints of P. This difference 
changes sign under reflection through P if b _L P and does not if b || P. 

Let be the cr-algebras of events that depend only on the portion of (rjb, K b )- 
configuration on Bp; explicitly J = a(rjb, b £ B^). Reflection positivity is, 
in its essence, a bound on the correlation between events (and random variables) 
from J ?p and ■'Pf- The precise definition is as follows: 


Definition 4.1. Let P be a probability measure on configurations {p b , Kb)bev, L and 
let E be the corresponding expectation. We say that P is reflection positive if for 
any plane of reflection P and any two bounded -measurable random variables 
X and Y the following inequalities hold: 


E (X9 P (Y)) = E (Y0 P (X)) (4.3) 


and 

E {XOp(X)) > 0. (4.4) 

Here, 0p(X) denotes the random variable X o Op. 
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Next we will discuss how reflection positivity underlines our principal tech¬ 
nical tool: chessboard estimates. Consider an event A that depends only on the 
( 77 b, Kb)-configurations on the plaquette with the lower-left corner at the torus ori¬ 
gin. We will call such an A a plaquette event. For each a; g Tj,,we define x) x (A) to 
be the event depending only on the configuration on the plaquette with the lower- 
left corner at x which is obtained from A as follows: If both components of x are 
even, then 0 X ( A) is simply the translate of A by x. In the remaining cases we first 
reflect A along the side(s) of the plaquette in the direction(s) where the component 
of x is odd, and then translate the resulting event appropriately. (Thus, there are 
four possible “versions” of t) x (A), depending on the parity of x.) 

Here is the desired consequence of reflection positivity: 

Theorem 4.2 (Chessboard estimate). Let P be a reflection-positive measure on 
configurations (r]bi Kb)beB L - Then for any plaquette events Ai ,..., A rn and any 
distinct sites x \,..., x m € Tp, 

m m 1 

p( n mao) < n p ( n itli • ^ 

3 =1 j =1 xgTi 

Proof. See [15, Theorem 2.2], 

The moral of this result—whose proof boils down to the Cauchy-Schwarz in¬ 
equality for the inner product X, Y 1 —> E(X0p(F))—is that the probability of any 
number of plaquette events factorizes, as a bound, into the product of probabilities. 
This is particularly useful for contour estimates (of course, provided that the word 
contour refers to a collection of plaquettes on each of which some “bad” event 
occurs). Indeed, by (4.5) the probability of a contour will be suppressed exponen¬ 
tially in the number of constituting plaquettes. 

In light of (4.5), our estimates will require good bounds on probabilities of the 
so called disseminated events DxeTi, Unfortunately, the event A is often 

a conglomerate of several, more elementary events which makes a direct estimate 
° f ru Ti •& X (A) complicated. Here the following subadditivity property will turn 
out to be useful. 

Lemma 4.3 (Subadditivity). Suppose that P is a reflection-positive measure and 
let Ai, A 2 , ■ ■ ■ and A be plaquette events such that A C [J • Aj. Then 

f( n ma ))*' < £ p( n («) 

zeT L j xST L 

Proof. This is Lemma 6.3 of [5]. 

Apart from the above reflections, which we will call direct , one estimate— 
namely (4.39)—in the proof of Theorem 2.2 requires the use of so called diagonal 
reflections. Assuming L is even, these are reflections in the planes P of sites of the 
form 


P± = {y € Tl: ei • (y-x) T e 2 • {y - x) G {0, L / 2 }}. 


(4.7) 
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Here a; is a site that the plane passes through and 81 and e > are the unit vectors 
in the x and (/-coordinate directions. As before, the plane has two components— 
one corresponding to 81 • (y — x) = ±82 • (y — x) and the other corresponding to 
81 • (y — x) = ±82 -(y — x) + l /2 —and it divides Tl into two equal parts. This puts 
us into the setting assumed in Definition 4.1. Some care is needed in the definition 
of reflected configurations: If b' is the bond obtained by reflecting b through P, 
then 


(°pv) b = 


Vb', 

-rib', 


if P = P+, 
if P = P_. 


(4.8) 


This is different compared to (4.1) because the reflection in P + preserves orienta¬ 
tions of the edges, while that in l b reverses them. 


Remark 4.4. While we will only apply these reflections in d = 2, we note that the 
generalization to higher dimensions is straightforward; just consider all planes as 
above with ( 81 , 82 ) replaced by various pairs ( 8 ,, 8 j) of distinct coordinate vectors. 
These reflections will of course preserve the orientations of all edges in directions 
distinct from 8 ,; and 8 j. 


4.2. Phase transitions on tori 

Here we will provide the proof of phase transition in the form stated in Theo¬ 
rem 2.4. We follow pretty much the standard approach to proofs of order-disorder 
transitions which dates all the way back to [9,20,21], A somewhat different ap¬ 
proach (motivated by another perspective) to this proof can be found in [7], 

In order to use the techniques decribed in the previous section, we have to de¬ 
termine when the extended gradient Gibbs measure Ql on obeys the conditions 
of reflection positivity. 

Proposition 4.5. Let V be of the form (2.5) with any probability measure g for 
which Z^y < 00 . Then Ql is reflection positive for both direct and diagonal 
reflections. 

Proof. The proof is the same for both types of reflections so we we proceed fairly 
generally. Pick a plane of reflection P. Let z be a site on P and let us reexpress 
the rib’s back in terms of the </>’s with the convention that (f> z = 0. Then 

^L(d? 7 b) = 6(d(j> z ) d <f> x . (4.9) 

X^Z 

Next, let us introduce the quantity 

W(r), k) = i K 'bVb + j KbT] b- ( 4 - 10 ) 

b€B +\P 

(We note in passing that the removal of P from the first sum is non-trivial even 
for diagonal reflections once d > 3.) Clearly, W is jFp -measurable and the full 
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( 77 , k )-interaction is simply W(rj, k) + ( 9pW)(r), k). The Gibbs measure Qp can 
then be written 

Q L (dr,,dn) = J— e - w ^M^ w )^M5(dc(, z ) ( J] d<^ [] p{d Kb ) 

Zl ' v Ket 

(4.11) 

Now pick a bounded, -measurable function X = X (//. k) and integrate the 
function X6pX with respect to the torus measure Qp. If is the a -algebra 
generated by random variables o x and n b , with x and b “on” P, we have 

E Ql (X0 p X\&p) oc (J X(V</>,K)e- w W' K ) J] df x J] e (d/c 6 )) > 0, 

' xGT +\P 6eB+\P / 

(4.12) 

where the values of (« b , <!> x ) on P are implicit in the integral. This proves the 
property in (4.4); the identity (4.3) follows by the reflection symmetry of Qp. 


Let us consider two good plaquette events, f/ or( j and <7 ( | IS , that all edges on the 
plaquette are ordered and disordered, respectively. Let B = (t) lird u t/ ( ii S ) c denote 
the corresponding bad event. Given a plaquette event A, let 


$L,p(A) 


q l { n m *)) 


(4.13) 


abbreviate the quantity on the right-hand side of (4.5) and define 


i{A) = limsup sup $l, p (A). (4.14) 

L—>00 0<p<l 


The calculations from Sect. 3 then permit us to draw the following conclusion: 


Lemma 4.6. For each 6 > 0 there exists c > 0 such that if Kq > ckd, then 


m < 5 . 


(4.15) 

Moreover, there exist po, p\ £ (0,1) such that 



limsup3L jP ((/ or d) < <5, 

L—>00 

P<Po, 

(4.16) 

and 

limsup3L, P (0dis) < S, 

L —»oo 

P > Pi- 

(4.17) 

Proof. The event B can be decomposed into a disjoint union of events B, each of 
which admits exactly one arrangement of ordered and disordered bonds around 
the plaquette; see Fig. 3 for the relevant patterns. If £>,; is an event of type a £ 
{O, D, UO, UD, MP, MA}, then 

limsup3z, )P (Bj) < exp{-[P L>a 

L —►oo 

(p) - min F & (p)]\. 
<5=0,D 

(4.18) 
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By Theorem 3.3, the right-hand side is bounded by 0(1)( k °/ Kd ) -1 / 8 , uniformly 
in p. Applying Lemma 4.3, we conclude that 3 l, p (B) is small uniformly in p £ 
[0,1] once L>1. (The values p = 0,1 are handled by a limiting argument.) 

The bounds (4.16-4.17) follow by the fact that 

F 0 (p) - F D (p) = - 21og—^- flog—, (4.19) 

1 —p Kd 

which is (large) negative for p near one and (large) positive for p near zero. □ 

From %{&) <C 1 we immediately infer that the bad events occur with very low 
frequency. Moreover, a standard argument shows that the two good events do not 
like to occur in the same configuration. An explicit form of this statement is as 
follows: 

Lemma 4.7. Let R°f d be the random variable from (2.12). There exists a constant 
C < oo such that for all (even) L > 1 and all p £ [0,1], 

E Ql {R°l\ 1 - K?j) < C iL , p (B). (4.20) 

Proof The claim follows from the fact that, for some constant C' < oo, 

Ql(MGo rd) n tySdis)) < C'iL, p (B)\ (4.21) 

uniformly in x. y £ T/. Indeed, the expectation in (4.20) is the average of the 
probabilities Ql(k 6 = Ko, = «d) over all 6, 6 £ M p . If x and y denotes 
the plaquettes containing the bonds b and b, respectively, then this probability is 
bounded by Q L (dx(G ord) n t?y(!?£ rd )). But C/^ rd = B U C/ d i s and so by (4.21) the 
latter probability is bounded by }l, p {B) + C'Il, p {B) 4 < (C' + 1 )$l,p(B), where 
we used $l, p (B) < 1 . 

It remains to prove (4.21). Consider the event t? x (l?ord) (T (^dis) where, with¬ 
out loss of generality, x y. We claim that on this event, the good plaquettes at x 
and y are separated from each other by a ^-connected circuit of bad plaquettes. To 
see this, consider the largest connected component of good plaquettes containing x 
and note that no plaquette neighboring on this component can be good, because 
(by definition) the events £/ or( j and () dls cannot occur at neighboring plaquettes (we 
are assuming that Ko ^ k d ). By chessboard estimates, the probability in Ql of any 
such (given) circuit is bounded by 3 l, p {B) to its size; a standard Peierls’ argument 
in toroidal geometry (cf the proof of [5, Lemma 3.2]) now shows that the proba¬ 
bility in (4.21) is dominated by the probability of the shortest possible contour— 
which is $l, p (B) 4 . (The contour argument requires that 3 l, p {B) be smaller than 
some constant, but this we may assume to be automatically satisfied because the 
left-hand side of (4.20) is less than one.) □ 

Now we are in a position to prove our claims concerning the torus state: 

Proof of Theorem 2.4. Let R°f d be the fraction of ordered bonds on T / (cf. (2.12)) 
and let Xl{p) be the expectation of R°f d in the extended torus state Ql with pa¬ 
rameter p. Since (1 — p) - l Bi lis log-convex in the variable h = log 
and 

01og((l -p)~^Z L y) 


\^l\xl(p) 


dh 


(4.22) 
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we can conclude that the function p i—> Xl{p) is non-decreasing. Moreover, as 
the thermodynamic limit of the torus free energy exists (cf Proposition 5.5 in 
Sect. 5.3), the limit xip) = liniL^oo Xl{p) exists at all but perhaps a countable 
number of p's —namely the set T> C [0,1] of points where the limiting free energy 
is not differentiable. 

Next we claim that QL{\R°i, d — xip) I > e o) tends to zero as L —» oo for 
all eo > 0 and all p (/-_ T>. Indeed, if this probability stays uniformly positive 
along some subsequence of L's for some eo > 0, then the boundedness of f?” d 
ensures that for some > 0 and some e > 0 we have Ql{R°I a > xip) + e) > C 
and Ql(R°[ a < xip) — e) > C f° r a il L i n this subsequence. Vaguely speaking, 
this implies p £ V because one is then able to extract two infinite-volume Gibbs 
states with distinct densities of ordered bonds. A formal proof goes as follows: 
Consider the cumulant generating function ^(/i) = IB/,1 -1 log EQ L (e h \ BL \ R L ) 
and note that its thermodynamic limit, <P(h ) = lim L^Kx>&L{h), is convex in h 
and differentiable at h = 0 whenever p (jL V. But Ql(R°[ a > xip) + e) > C in 
conjunction with the exponential Chebyshev inequality implies 

d> L (h) - h( X {p) + e) > (4.23) 


which by taking L —> oo and h , | 0 yields a lower bound on the right derivative 
at origin, > x(p) + e - By the same token Ql(R°l A < x(p) — e) > ( 

implies an upper bound on the left derivative, < x(p) ~ £• Hence, both 

probabilities can be uniformly positive only if p £ V. 

To prove the desired claim it remains to show that x jumps from values near 
zero to values near one at some p t € (0,1). To this end we first observe that 

lim E Ql (i?f (1 - R°l A )) = X (P) [1 - X(P)], P t V. (4.24) 

J-j —^OO 

This follows by the fact that on the event {xip) — e < < xip) + e l—whose 

probability tends to one as L —> oo—the quantity f?™ d ( 1 — f?™ d ) is bounded 
between [x(p) + e] (1 — xip) + e ) an d [xip) — e ] (1 — xip) ~ e ) provided e < 
min{x(p), 1 — xip)}- Lemma 4.7 now implies 

Xip) [1 - Xip)} < (4.25) 

with 3 (B) defined in (4.14). By Lemma 4.6, for each 6 > 0 there is a constant c > 0 
such that 

x(p)e [0,<5] U [1-5,1], P^V, (4.26) 

once k °/ Kd > c. But the bounds (4.16—4.17) ensure that x(p) G [0, <5] for p <C 1 
and\(p) £ [1 — (5,1] for 1— p <C 1. Hence, by the monotonicity of p i—> ^(p),there 
exists a unique value p t £ ( 0 , 1 ) such that xip) < 8 for p < p t while xip) >1 — 5 
for p > p t . In light of our previous reasoning, this proves the bounds (2.13-2.14). 
□ 
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4.3. Phase coexistence in infinite volume 


In order to prove Theorem 2.2, we will need to derive a concentration bound on 
the tilt of the toms states. This is the content of the following lemma: 


Lemma 4.8. Let A C Tl and let B/i be the set of bonds with both ends in A. Given 
a configuration (rib)beK L , we use Ua = Ua(v) to denote the vector 



6eB^ r nB^ 



£ ») 


(4.27) 


of empirical tilt of the configuration r/t, in A. Suppose that K mm = inf supp q > 0. 
Then 

Pl(\U a \> 5 ) < 4e _ 5 Kmin<52 l B/1 l (4.28) 

for each 6 > 0, each A C Tl and each L. 


Proof. We will derive a bound on the exponential moment of IJ,\. Let us fix a 
collection of numbers (vb)beM L € and let Ql,(ki ,) be the conditional law of 

the rf s given a configuration of the n's. Let Ql . m in be the corresponding law when 
all Kb = K m i n • In view of the fact that Ql,(k,,) and Ql .min are Gaussian measures 
and Kb > h' m in, we have 

Var ^,(* 6 ) 

v bGB L ' v bG B l 


( Yl Vbllb ) - Var Oi,»in ( Yl Vbllb ) • (4 - 29) 


(Note that both measures enforce the same loop conditions.) The right-hand side 
is best calculated in terms of the gradients. The result is 


Var <^,W) ( 


Y VbT l b 

bGBi 


< 


1 

b'min 


bGB L 


vl- 


(4.30) 


The fact that Eq l (ijb) — 0 and the identity £(e x ) = e Bj,(:+ 2 Var ( x ), valid for 
any Gaussian random variable, now allow us to conclude 


E Q L ( exp{ Y v bPb} ) < ex p{^- Y v l}- 

A bGB L ' mm bGB l 


(4.31) 


Choosing Vb = A • 


on B/| and zero otherwise, we get 


(4.32) 


Noting that \Ua\ > 5 implies that at least one of the components of Ua is larger 
(in absolute value) than the desired bound follows by a standard exponetial- 
Chebyshev estimate. 
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Remark 4.9. We note that the symmetry of the law of the 77’s in Ql.(k 6 ) is crucial 
for the above argument. In particular, it is not clear how to control the tightness of 
the empirical tilt Ua in the measure obtained by normalizing exp{— Yhb V(Vb + 
^6)}t / i(dt?), where hb = h ■ b is a “built-in” tilt. In the strictly convex cases, these 
measures were used by Funaki and Spohn [17] to construct an infinite-volume 
gradient Gibbs state with a given value of the tilt. 


Proof of Theorem 2.2. With Theorem 2.4 at our disposal, the argument is fairly 
straightforward. Consider a weak (subsequential) limit of the torus states at p > p t 
and then consider another weak limit of these states as p J, p t . Denote the result 
by /2 0r d. Next let us perform a similar limit as p | p t and let us denote the resulting 
measure by /7dis- As is easy to check, both measures are extended gradient Gibbs 
measures at parameter p t . 

Next we will show that the two measures are distinct measures of zero tilt. 
To this end we recall that, by (2.14) and the invariance of Ql under rotations, 
liminfi^oo = k 0 ) > 1 — e when p > p t while (2.13) implies that 

limsup L _ >0O = kq) < e when p < p,. But {/p, = k 0 } is a local event 

and so 

Ford( k& = «o) > 1 - e (4.33) 

while 

Fdis (Kb = Ko) < e, (4-34) 

for all b; i.e., /i or d pidis- Moreover, the bound (4.28)—being uniform in p and L — 
survives the above limits unscathed and so the tilt is exponentially tight in volume 
for both /2o r d and /idi s . It follows that Ua —> 0 as A | Z 2 almost surely with respect 
to both /i or d and /T d is; he., both measures are supported entirely on configurations 
with zero tilt. 

It remains to prove the inequalities (2.10-2.11) and thereby ensure that the //- 
marginals p 0 rd and pdis of Ford and /Idis, respectively, are distinct as claimed in the 
statement of the theorem. The first bound is a consequence of the identity 

lim E QL (K b r]l) = (4.35) 

L—> 00 4 


which extends via the aforementioned limits to p or d (as well as Fdis)- Indeed, using 
Chebyshev’s inequality and the fact that Ford(^b = «d) < e we get 


Ford(bt’0Pfr ' ) r Ford(tvoPb ^ 5 tvf, — h'.()j 


< Ford (KbPb > A 2 ) < 


E UtMbvl) = J_ 

A 2 4A 2 ■ 


To prove (4.35), the translation and rotation invariance of () /, gives us 
e Ql (KbFb) = E Ql (eQlm ( |^| Yl Kb 7 lb) ) • 


(4.36) 


(4.37) 


Let Zl denote the integral of exp{—4 Y^b K bVb } with respect to fl- Since we 

have rj) = /3l Ti ^ 1 ^i(d?7), simple scaling of all fields yields Z L ^p Kb ^ = 
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j3 2 (l T r*I ^Zl,( K b )- Intepreting the inner expectation above as the (negative) / 3- 
derivative of at (3 = 1, we get 


Eq 





ITlI-i 

2|B l | 


(4.38) 


From here (4.35) follows by taking L —> oo on the right-hand side. 

As to the inequality (2.11) for the disordered state, here we first use that the 
diagonal reflection allows us to disseminate the event {KbVb — ^ 2 } around any 
plaquette containing b. Explicitly, if ( 6 i, 62 , 63 , 64 ) is a plaquette, then 

QLiKb^ < A 2 ) < Q l ( P| {K b ril < \ 2 }^j . (4.39) 


(We are using that the event in question is even in // and so the changes of sign 
of rjb are immaterial.) Direct reflections now permit us to disseminate the resulting 
plaquette event all over the torus: 

< A 2 ) < Ql( P {K~ b vl < A 2 }) . (4.40) 

be s L 

Bounding the indicator of the giant intersection by 

e^-D^IB.I exp{-i(/3 - 1) £ W7&}, (4.41) 


for /? > 1 , and invoking the scaling of the partition function Z L ^ Kb y we deduce 


QL^bVb < A 2 ) < 


A/3-l)A 2 |Bi| 


P 


i(|1Tz,|-l) 




(4.42) 


Choosing (3 = A 2 , letting L —> 00 and p | p t , we thus conclude 

Pdi s (nbVb < A 2 ) < C 1 A 14 . (4.43) 

Noting that /idis(ft&? 7 b — ^ 2 ) — AGis^D ?? 2 — ^ 2 ) — e - the bound (2.11) is also 
proved. □ 


5. Duality arguments 

The goal of this section is to prove Theorem 2.5. For that we will establish an 
interesting duality that relates the model with parameter p to the same model with 
parameter 1 — p. 
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5.1. Preliminary considerations 


The duality relation that our model (1.6) satisfies boils down, more or less, to 
an algebraic fact that the plaquette condition (1.3), represented by the delta func¬ 
tion S{rj bl + r] b2 — r]b 3 — can formally be written as 

S(v bl +Vb 2 - Vb 3 - = J ^e‘^^1+^-^3-^). (5.1) 

We interpret the variable <f>* as the dual field that is associated with the plaquette 
{y]b\ i Vb 2 i Vb 3 , Vb 4 )- As it turns out (see Theorem 5.3), by integrating the rj's with 
the cf>*'s fixed a gradient measure is produced whose interaction is the same as for 
the ry’s, except that the n b 's get replaced by l/n b ’s. This means that if we assume 
that 

k 0 k d = 1, (5.2) 

which is permissible in light of the remarks at the beginning of Section 2.2, then the 
duality simply exchanges n 0 and re D ! We will assume that (5.2) holds throughout 
this entire section. 

The aforementioned transformation works nicely for the plaquette conditions 
which guarantee that the rj's can locally be integrated back to the f’s. However, 
in two-dimensional torus geometry, two additional global constraints are also re¬ 
quired to ensure the global correspondence between the gradients rj and the f’s. 
These constraints, which are by definition built into the a priori measure vl from 
Sect. 2, do not transform as nicely as the local plaquette conditions. To capture 
these subtleties, we will now define another a priori measure that differs from vl 
in that it disregards these global constraints. 

Consider the linear subspace Xf D Xl of R Bi that is characterized by the 
equations ij bl +rj b2 — fjb 3 — % 4 = 0 for each plaquette (bi,b 2 , b 3 , 64 ). This space in¬ 
herits the Euclidean metric from R Bi ; we define vf as the corresponding Lebesgue 
measure on Xf scaled by a constant Cl which will be determined momentarily. 
In order to make the link with vl , we define 


r/verl — ^ ' Vx+ei and TJhor — y ' TJx +e 2 • (5.3) 

xCLTl xETl 

Clearly, 

Xl = {v € XI: ?? vert = 0,77hor = 0}. (5.4) 

Consider also the projection IIi.: Xf —> Xl which is defined, for any configura¬ 
tion p G R Bl , by 


{n L v)b 


Tjb ^/vert 5 

?]b "JJI ^Thor 1 


if b € 

if fee B^ or . 


(5.5) 


Then we have: 
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Lemma 5.1. There exist constants Cl such that, in the sense of distributions, 
v*l{ d?y) = (l 2 f dd 5(j] bl + ?/ &2 - r/ b3 - p bi - 9) d? ?& . (5.6) 

' 1/111 ( 61 , 62 , 63 , 64 ) ' 6 eTi 

Moreover, we have 

U L (dr]) = 6'i(d77)5(77hor)()(?7vert) (5.7) 

and 

isl(dri) = v L o i 7 L (d? 7 ) An^dr?). (5.8) 

//ere, (dr/) is a multiple of the Lebesgue measure on the two-dimensional 
space Tlf (0) n which can be formally identified with d?/hoi-d?/ V ert- 

Proof We begin with (5.6). Consider the orthogonal decomposition RI Bl I = Xf ® 
(Xf) ± . Clearly, dim Xf = L 2 + 1. Choosing an orthonormal basis w\ ,..., w n 
in (Xf)- 1 (where n = dim(A’£)- L = L 2 — 1) the measure v* L can be written as 

n 

o* L (dp) =cl( n • vj) n dr ? 6 • (5 - 9) 

j=i 6 eBi 

Let £ n denote the vectors in Rl Bi l such that if 7 r = (b \, 62 , 63 , 64 ) then ■ p = 
r] bl + Ph 2 ~ V b 3 — Vbi- Then £ T £ {Xl) ± with all but one of these vectors linearly 
independent. This means that we can replace the linear functionals // i—^Wj-rj by 
the plaquette conditions. Fixing a particular plaquette, 7 To, we find that 

"so dr ?)=(n ■ n) ] n &r ib (5.io) 

^■71^710 6gT_l 

provided that 

Cl = |det(tUj • £^)\ = \Jdet(£ n ■ £ n >). (5.11) 

The expression (5.10) is now easily checked to be equivalent to (5.6): Applying 
the constraints from the plaquettes distinct from 7r, we find that rj bl + rj b2 — rj b3 — 
r] bi = (1 — L 2 )9. The corresponding (5-function becomes S(L 2 9), and so we can 
set 9 = 0 in the remaining (5-functions. Integration over 9 yields an overall multi¬ 
plier f R 5(L 2 9)d9 = \IL 2 . 

In order to prove (5.7), pick a subtree T of T l as follows: T contains the 
horizontal bonds in {£»i + £e\: £ = 0,..., L — 2} and the vertical bonds in 
{62 + £e\ + me 2 : £,m = 0,..., L — 2}. As is easy to check, T is a span¬ 
ning tree. Denoting by 7 l the measure on the right-hand side of (5.7) pick a 
bounded, continuous function /: Rl Bi —> R with bounded support and consider 
the integral f / ( 77 ) 74 , (d?/). The complement of T contains exactly L 2 + 1 edges 
and there are as many (5-functions in (5.10) and (5.7), in which all rj b , b ^ T, 
appear with coefficient ±1. We may thus resolve these constraints and substi¬ 
tute for all {rj b : b ^ T} into /—call the result of this substitution f{rj). Then 
we can integrate all of these variables which reduces our attention to the inte- 
gral / f(j]) n facT drj b . 
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As is easy to check, the transformation r/i, = <p y — (f> x for b = ( x, y) with 
the convention <f >o = 0 turns the measure Tiber ^Vb into <5(d<^o) Wx^lq^x and 
makes f(rj) into /(V(/>). We have thus deduced 

/ /(r?)7r(dr?) = [ /(r?) d?? & = / /(V0)<5(d^ o ) J|d^ x . 

Jr\t\ ££ J r \t l \ 

(5.12) 

From here we get (5.7) by noting that the latter integral can also be written as 

J7(>7H(d»/)- 

To derive = vl ° IILdi% 0 idr] WT t, let us write Xf = Xl © {II L 1 (0) (T Xf). 
Since is the C^-multiple of the Lebesgue measure on Xj and since 77h or and r/ velt 
represent orthogonal coordinates in 11^(0) D X£, we have 

v* L {drj) = C l \ Xl ° n L (dr)) L^d^ L~ 2 dri vert , (5.13) 

where Xx L is the Lebesgue measure on Xl- Plugging into (5.7) we find that vl = 
ClL~ 4 Xx l which in turn implies (5.8). □ 

Remark 5.2. It is of some interest to note that the measure z/£ is also reflection pos¬ 
itive for direct reflections. One proof of this fact goes by replacing the (5-functions 
in (5.6) by Gaussian kernels and noting that the linear term in 6 (in the exponent) 
exactly cancels. The status of reflection positivity for the diagonal reflections is 
unclear. 


5.2. Duality for inhomogeneous Gaussian measures 


Now we can state the principal duality relation. For that let T£ denote the dual 
torus which is simply a copy of shifted by half lattice spacing in each direction. 
Let B>T denote the set of dual edges. We will adopt the convention that if b is a 
direct edge, then its dual—i.e., the unique edge in B£ that cuts through b —will be 
denoted by b *. Then we have: 


Theorem 5.3. Given two collections (bCb)beB L and (bt^beBt of positive weights 
on Ml, consider the partition functions 

Z L,(K b ) = / vl{ dr]) expj —- ^ btb?7 2 } (5.14) 


and 

z l,( k *) = J u t(dv) exp 

If (Kb)bES L an d (Kj)beB L are dual in the sense that 


{"2 ^ 4%2 }' 


bPBr, 


then 


Kb 


b £ Ml, 


z U*t) = 


n 


bSBi 


j L,(k b )- 


(5.15) 


(5.16) 

(5.17) 
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Proof. We will cast the partition function Z) (Kt . f into the form on the right-hand 
side of (5.17). Let us regard this partition function as defined on the dual torus T£. 
The proof commences by rewriting the definition (5.6) with the help of (5.1) as 


H( d ?7) = L 


dO / [] 

; J ^xJ L 


2i r 


exp {iE 4>x(Tln*(x) ~ #)} J II d 7b*. 

x€T l ' 6 * 68 ^ 

(5.18) 


where T) n */ X ) is the plaquette curl for the dual plaquette 7r* (x) with the center at x. 
Rearranging terms and multiplying by the exponential (Gaussian) weight from 
(5.15), we are thus supposed to integrate the function 


d9 


n 

'x£Yl 


d <p x 

27r 


exp 


He 




b*eM* 


- 2ir]i,*\7 b (f>) - W 

x£T l 

(5.19) 


against the (unconstrained) Lebesgue measure d%. Here Nb4> = <j) y — (j) x 

if b = (x , y) is dual to the bond b*. Completing the squares and integrating over 
the p’s produces the function 


L 2 


n 

b*e B* 


b* -I 


d6 


ImJl 


P dfx 

xeT L 


exp 


--vi 

2 b{tt L K b* 


(v b </>) 2 -w V </> : 


xeT L 

(5.20) 


Invoking (5.16), we can replace all 1 /nf by n b . The integral over 6 then yields 27 t 
times the ^-function of <f> x which—by the substitution <j> x i—> (f> x + ytjj&o that 

has no effect on the rest of the integral—can be converted to 2tt6(<Po)- Invoking 
the definition (2.1) of vp- this leads to the partition function (5.14). □ 


Remark 5.4. Let Ql. p be the extended gradient Gibbs measure Ql for g = pS Ko + 
(1 — p)6 KD with parameter p and let Q* L be the corresponding measure with the 
a priori measure vp replaced by i/* L . Then the above duality shows that the law 
of (Kb) governed by Ql.p is the same as the law of its dual (k£) —defined via 
(5.16)—in measure Q* L , once p and p* are related by 


P P* 

1 —p 1 — p+ 



(5.21) 


Indeed, the probability in measure Q* L of seeing the configuration K) with N£ 
ordered bonds and TVq disordered bonds is proportional to pH (1 — p*) N ° Z* L y 
Considering the dual configuration («(,) and letting No = Nq denote the number 
of disordered bonds and No = N£ the number of ordered bonds in (Kb), we thus 
have 


pH(l - Pit )^Zl A<) = 27rL 2 (p*^) JVD (( 1 -p*)V^) N °Z L , {Kb y (5.22) 

For p and p* related as in (5.21), the right-hand side is proportional to the prob¬ 
ability of (Kb) in measure Ql.p- We believe that the difference between the two 
measures disappears in the limit L —» oo and so the /v-marginals of the states /j or< i 
and /idi S at p t can be considered to be dual to each other. However, we will not 
pursue this detail at any level of rigor. 
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5.3. Computing the transition point 

In order to use effectively the duality relation from Theorem 5.3, we have to show 
that the difference in the a priori measure can be neglected. We will do this by 
showing that both partition functions lead to the same free energy. This is some¬ 
what subtle due to the presence (and absence) of various constraints, so we will 
carry out the proof in detail. 

Proposition 5.5. Let ql( dft) = J~[ beBt \P^o(^ K b) + (1 — p)£ KD (dK&)] and recall 
that Zl,v denotes the integral of with respect to ql- Similarly, let Z* L v 

denote the integral of Z* L with respect to ql- Then (the following limits exist 
as L —> oo and) 

lim \ogZ Ly = lim E-j-log Z* L V (5.23) 

L-> oo | 1 l| L—>oo |I l | 

for all p £ [ 0 , 1 ]. 

Before we commence with the proof, let us establish the following variance 
bounds for homogeneous Gaussian measures relative to the a priori measure vj. 
and v^'. 

Lemma 5.6. Let /// he the (standard) Gaussian gradient measure 

Pl( dt?) oc exp|-i Vb}vL(dri) (5.24) 

be E l 

and p,* L be the measure obtained by replacing Ul by uf. For m = 1,..., L, let 

m —1 

= ^ ^ ^7(^ei,(^+l)ei) • (5.25) 

i=0 

There exists an absolute constant C 3 > 0 such that for all L > 1 and all m = 
1 

Var Mi (Y m ) < Var M * (Y m ) < c 3 (l + logm). (5.26) 


Proof. In measure /j/\, we can reintroduce back the fields (4> x ) and Y m then equals 
(/) m gi. Discrete Fourier transform implies that 


Var u. L {(j) x ) 


1 

If 


E 

kei L 


|1 - e ix ' fc | 2 

D(k) 


(5.27) 


where is the reciprocal torus and D(k) = |1 — e lfcl | 2 + |1 — e lfc2 | 2 is the 
discrete (torus) Laplacian. Simple estimates show that the sum is bounded by a 
constant times 1 + log \x\, uniformly in L. Hence, Var^ ( Y m ) < c 3 (l + log?n) 
for some absolute constant c 3 . 

As for the other measure, we recall the definitions (5.3) and use these to 
write ??b = Vb</> + -p-??hor if b is horizontal (and i% = Vb</> + -p-?? V ert if b is verti¬ 
cal). The fact that the Gaussian field is homogeneous implies—via (5.8)—that the 
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fields (<f> x ) and the variables rj ven and t^or are independent with ((f> x ) distributed 
according to p,L and r] ven and ? 7 h or Gaussian with mean zero and variance L ' 2 / 2 - In 
this case Y m = 0 mgl + and so we get 

777,2 

Var^t* (Y m ) = Var ML (Y m ) + (5.28) 

But m < L and so the correction is bounded for all L. □ 

Proof of Proposition 5.5. The proof follows the expected line: To compensate for 
the lack of obvious subadditivity of the torus partition function, we will first relate 
the periodic boundary condition to a “fixed” boundary condition. Then we will 
establish subadditivity—and hence the existence of the free energy—for the latter 
boundary condition. 

Fix M > 0 and consider the partition function 44 defined as follows. Let Al 
be a box of Lx L sites and consider the set of edges with both ends in Al. Let 
Vf M ' > (dr)) be as in (2.1) subject to the restriction that \<p x \ < M for all x on the 
internal boundary of Al - Let 



We will now provide upper and lower bounds between the partition functions Z^y 
(resp. Zf v ) and Zfll , for a well defined range of values of M. 

Comparing explicit expressions for Zl.v and 44 and using Kb < Ko, we 
get 

Z L ,v > 4^ exp{-i Ko (2M) 2 (2 L)}. (5.30) 

To derive an opposite inequality, note that for Kb > kd we get that VarQ L {k ) (<f> x ) > 
Var^^ (4>x)/kd, where is as in (5.24). Invoking one more time the Gaussian 
identity E(e x ) = e BJf + 5 Var (- Y ) in conjunction with Lemma 5.6, yields 

QM, > M) < exp {-i—(5.31) 

Hence, if M log L we have that with probability at least V 2 in measure Ql, all 
variables <f> x are in the interval [— M, M]. Since the interaction that wraps Al into 
the torus is of definite sign, it follows that 


Z L ,v < 2 44 


(5.32) 


for all L and all M 5P log L. 

Concerning the star-partition function. Lemma 5.6 makes the proof of (5.32) 
exactly the same. As for the alternative of (5.30), we invoke (5.8) and restrict 
all | (j> x | on the internal boundary of Al to values less than M and |? 7 hor| and r/ vcn 
to values less than ML 2 . Since \r]b\ = \Vb<t> + j?Vvcn < 2 M + M = 3 M for 
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every vertical bond that wraps Al into the toms (and similarly for the horizontal 
bonds), we now get 

Zi tV > Z^y (2ML) 2 exp{ — 5 k 0 ( 3 M) 2 ( 2 L)}, (5.33) 

where the factor (2 AIL) 2 comes from the integration over rj vea and j^or- We con¬ 
clude that, for log L M = oip/L), the partition functions Zl.v, Z* L v and zj^y 
lead to the same free energy, provided at least one of these exists. 

It remains to establish that the partition function z\ M y is (approximately) sub- 
multiplicative for some choice of M = Ml. Choose, e.g.. Ml = (logL ) 2 and 
let p > 1 be an integer. If two neighbors have their <j)’s between —Ml and Ml, 
the energy across the bond is at most ^Ko(4Ml) 2 ■ Splitting A p l into p 2 boxes of 
size L, and restricting the </>’s to [—Ml, Ml] on the internal boundaries of these 
boxes, we thus get 

> [Z ( l M v l Y exp{-i Ko (2 M l ) 2 2(p - 1 )L}. (5.34) 

The exponent can be bounded below by ( pL) 3/2 — p 2 L 3/2 = — {p 2 — p 3/2 )L 3/2 
for L sufficiently large which implies that p i—> [Z^ p y' > exp { — {pL) 3,2 }] 1 ^ pL ^ 
is increasing for all p > 1 and all L 1. This proves the claim for limits along 
multiples of any fixed L; to get the values “in-between” we just need to realize 
that, as before, Z^+^y 1 > Z^y^ e°( kLM l\ for any fixed k. □ 

Now we finally prove our claim concerning the value of the transitional p: 

(p) 

Proof of Theorem 2.5. Let Z ]v denote the integral of Zl ,( Kb ) with respect to the 

a priori measure ql{ dtc) with parameter p and let Z^fy denote the analogous 
quantity for Z* L The arguments leading up to (5.22) then yield 

Zffi = Z^(27rL 2 )(p^+( 1 -p*)^)' Bil (5-35) 

whenever p* is dual to p in the sense of (5.21). Thus, using F(p) to denote the 
limit in (5.23) with the negative sign, we have 

F(p*) = F(p) - 21og(p* v /Kjj + (1 -p+)y/KE). (5.36) 

Now, as a glance at the proof of Theorem 2.4 reveals, the value p t is defined as the 
unique point where the derivative of F(p), which at the continuity points of p i—> 
x(p) is simply F'(p) = 2 xip) — 1, jumps from values near —1 to values near +1. 
Eq. (5.36) then forces the jump to occur at the self-dual point p t = p. In light 
of (5.21), this proves (2.15). □ 
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